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Compaction  with  Automatic  Jog 
Introduction 

F.  Milter  Matey 

Laboratory  for  Computer  Science 
Massachusetts  Institute  of  Technology,  Cambridge 


Abstract— This  paper  presents  a  novel  polynomial* time  algorithm  for  com* 
parting  a  VLSI  layout.  Compared  to  previous  algorithms,  the  algorithm 
promises  to  produce  higher  quality  output  while  reducing  the  need  Cor  de¬ 
signer  intervention.  The  performance  gain  is  realized  by  converting  wires 
into  constraints  on  the  positions  of  the  active  devices.  These  constraints  can 
be  solved  by  graph-theoretic  techniques  to  yield  optimal  positions  for  chip 
components.  A  single-layer  router  is  then  used  to  restore  the  wires  to  the 
layout,  using  as  many  jogs  as  necessary. 

1.  Introduction 

An  automated  compaction  procedure  is  an  effective  tool  for  cutting  produc¬ 
tion  costs  of  a  VLSI  circuit  at  low  cost  to  the  designer,  because  the  yield  of 
fabricated  chips  is  strongly  dependent  cm  the  total  circuit  area-^An  effective 
compaction  system  also  reduces  design  time  by  freeing  the  designer  from 
continual  concern  over  design  rules.  If  excess  layout  space  can  be  removed 
automatically,  the  designer  can  sketch  a  layout  without  making  constant  ef¬ 
forts  to  conserve  area.  For  these  reasons,  compaction  algorithms  have  gained 
widespread  attention  in  the  VLSI  literature  [4,  5,  9,  11],  and  have  been  in¬ 
corporated  into  many  recent  computer-aided  circuit  design  systems  [2, 4, 10, 


Most  compaction  algorithms,  including  the  one  described  here,  com¬ 
press  a  layout  in  one  dimension  only.  To  reduce  both  dimensions,  the  layout 
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can  be  alternately  compacted  in  x  and  y  until  the  algorithm  can  find  no  fur¬ 
ther  improvement.  Compaction  in  two  dimensions  simultaneously  is  theoret¬ 
ically  difficult,  although  it  seems  to  work  well  in  practice  [5].  The  remainder 
of  this  paper  assumes  that  the  direction  of  compaction  is  horisantal. 

1.1.  Constraint-butd  compaction 

The  more  recent  systems  |4,10]  use  a  conttroint-bcucd  technique  to 
provide  one-dimensional  compaction.  These  procedures  begin  by  assigning 
to  each  layout  component  i  a  variable  x,,  which  represents  the  x-coordinate  of 
that  component  s  leftmost  point.  The  detign  rules  at  the  fabrication  process 
are  then  used  to  derive  constraints  on  the  positions  of  the  components.  For 
example,  if  device  i  lies  to  the  left  of  device  and  such  devices  must  remain 
at  least  2  units  apart  in  order  to  function  reliably,  the  compactor  generates 
a  constraint  xy  —  x,  >  2  +  to,,  where  u>,  is  the  width  of  component  t.  (We 
make  the  usual  assumption  that  components  are  not  allowed  to  jump  over 
one  another.) 

The  design  rules  lead  naturally  to  a  set  of  constraints  with  very  nice 
properties.  First  of  all,  the  constraints  are  not  especially  difficult  to  compute 
[9].  Second,  they  are  sufficient  to  guarantee  that  the  compacted  layout  is 
legal.  Third,  all  constraints  are  timplt  Untar  incgualitiei:  they  have  the 
form 

*i  -  *i  >  «v  , 

where  x,  and  xy  are  two  of  the  variables  assigned  to  layout  components,  and 
o,y  is  a  constant. 

Because  of  the  simple  form  of  the  inequalities,  they  can  be  solved 
efficiently  by  graph-theoretic  techniques.  One  constructs  an  edge-weighted 
graph  whose  vertices  represent  the  variables  x,,  and  in  which  an  edge  of 
weight  a,y  from  node  x,  to  node  xy  represents  the  constraint  xy  -  x<  >  o,y. 
An  assignment  to  the  variables  x<  that  satisfies  all  the  constraints  is  then 
determined  by  a  longest- path  computation  on  the  graph.  The  resulting 
values  specify  the  optimal  positions  of  the  components  in  the  compacted 
layout.  A  good  introduction  to  constraint-based  compaction  may  be  found 
in  [5];  common  algorithms  for  computing  longest  paths  are  discussed  in  [8]. 

1.2.  Automatic  jog  introduction 

In  order  to  perform  any  sort  of  compaction,  the  components  of  the 
layout  must  be  differentiated  into  modules,  which  are  fixed  in  sise  and  shape, 
and  wires,  which  are  flexible.  Common  procedures  for  generating  design 
rule  constraints  (4,5,9]  assume  that  wires  are  simply  rectangular  regions  of 
variable  length,  and  otherwise  identical  to  modules.  A  vertical  wire,  for 
example,  would  be  assigned  an  x-coordinate  during  horisontal  compaction, 


and  could  only  be  moved  rigidly  from  aide  to  ride.  But  one  would  often  ike 
*  previously  straight  wire  to  bend  mound  an  obstacle  during  compaction,  if 
the  area  at  the  circuit  could  thereby  be  reduced. 

This  problem  is  not  eerily  overcome.  Many  eye terns  [4,  17]  attempt 
to  solve  it  by  allowing  the  deeigurr  to  speedy  joy  paints  at  which  wises 
may  bend.  Compaction  then  becomes  an  interactive  procedure  in  which  the 
designer  repeatedly  examinee  the  compacted  layout,  adds  mace  potential  Jog 
points,  and  retries  the  compaction  operation.  Other  systems  (4]  attempt  to 
insert  jogs  automatically,  using  ad  hoe  techniques  which  are  not  gnamnteod 
to  be  effective.  One  technique  that  will  work  is  to  insert  aU  possible  jag  paints 
into  every  wire.  Tbe  resulting  algorithm  baa  the  theoretical  drawback  that 
it  requires  exponential  time  on  some  inputs;  even  in  practice,  it  is  fikaiy  to 
consume  large  amounts  of  time  and  memory. 

The  polynomial-time  algorithm  presented  in  this  paper  has  the  capa¬ 
bility  to  introduce  jog  points  automatically  m  an  optimal  fashion.  It  can  thus 
be  expected  to  produce  high  quality  output  with  little  designer  intervention 
Automatic  jog  introduction  is  achieved  by  treating  wine  not  as  solid  objects, 
but  only  as  indicators  of  the  topology  of  the  layout.  Constraints  between 
modules  no  longer  express  design  rules  directly;  instead,  they  will  ensure 
that  there  exist  paths  for  the  wires,  having  the  given  topology,  that  sat¬ 
isfy  the  design  rules.  The  new  constraints,  caOed  rentability  conditions,  can 
be  formulated  as  simple  linear  inequalities,  and  solved  as  usual.  When  the 
optimal  module  placements  have  been  established,  the  new  wire  paths  are 
determined  by  a  ringJe- layer  router  presented  in  |6].  This  router  generates 
no  "unnecessary  U’s,"  and  therefore  minimises  wire  lengths,  given  that  the 
layout  topology  is  fixed. 


This  approach  to  compaction  depends  m  >he  ability  to  | 
plete  routability  conditions  for  a  layout.  Until  recently,  each  conditions i 
known  only  for  curtain  channel  routing  problems  [7,16].  The  present  work  is 
made  paaaible  by  tbe  theory  of  planar  mutability  developed  in  [1]  and  [14]. 


Sr.( 


I  I S  CL+\  i  a  *  4  f,  -i 


J 


1.3.  Organisation  o/  tkt  paper 

Tbe  remainder  of  this  paper  is  organised  as  follows.  I>ectien2  states  the 
definitions  and  theoretical  results  that  underlie  the  new  compaction  method. 
In  this  extended  abstract,  ell  proofs  ere  emitted  Sections  show*  bow  the 
circuit  layout  is  converted  to  a  data  structure  appropriate  for  compaction, 
and  Section  4  details  tbe  body  of  tbe  compaction  algorithm 
Section  Vieveral  improvements  to  the  algorithm,  that  should  make  H  run 
considerably  faster.  Section  6  comments  briafiy^qn  the  algorithm’s  correct¬ 
ness  proof.  Finally;  I  conclude  in?  Secticrfc  some  extensions  af-my> 

results,  and  a  discussion  of  the  practical  value  of  et^com paction  algorithm. 
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2.  Sketches  and  planar  routability 

The  principal  data  atructurc  used  by  the  compaction  algorithm  ia  called  a 
sketch.  A  sketch  represents  one  layer  of  a  VLSI  circuit,  including  both  Axed 
object!  and  the  topology  of  the  interconnecting  wirea.  The  algorithms  in  this 
paper  process  only  one  sketch  at  a  time,  without  loaa  of  generality.  Com* 
paction  of  a  circuit  with  multiple  layer*  can  be  accomplished  by  computing 
the  constraints  far  each  layer  independently,  and  merging  the  constraint  sys¬ 
tems.  This  section  defines  precisely  what  1  mean  by  a  sketch,  and  states  the 
theorem  from  [6]  that  determines  routability  conditions  for  a  sketch.  For 
simplicity,  the  model  we  will  use  is  more  restricted  than  necessary.  Section  7 
notes  that  some  of  the  restrictions  can  be  relaxed. 

2.1.  Definition  o/  a  sketch 

A  sketch  is  an  ordered  pair  (F,W)  consisting  of  a  finite  set  F  of  /* a- 
(ures,  which  are  prints  and  straight  line  segments,  and  a  finite  set  IV  of 
wires,  which  are  simple  paths  in  the  plane.  Figure  1  allows  an  of  a 

■ketch.  Modules  are  represented  as  collections  of  features,  because  far  tech¬ 
nical  reasons,  terminals  must  be  separated  from  other  features.  The  features 
and  wires  of  a  sketch  must  satisfy  the  following  restrictions: 

(1)  Distinct  components  of  the  sketch  (feature*  and  wires)  may  intersect 
only  at  their  endpoints. 

(2)  Distinct  wires  may  not  intersect,  and  no  wire  may  cross  itself. 

(3)  Each  wire  touches  exactly  two  features,  which  are  single  prints  lying 
at  the  endpoints  of  the  wire.  They  are  called  the  terminals  of  the  wire. 
(Multitenninal  wires  can  also  be  handled;  see  Section  7.) 

(4)  Four  of  the  features  of  the  sketch  farm  a  bounding  box  around  the  other 
components. 

When  referring  to  “points  in  the  sketch,"  we  will  mean  points  lying  on  fea¬ 
tures  in  the  sketch.  The  connected  groups  of  features  are  called  the  obstacle# 
of  the  sketch.  The  definitions  imply  that  each  terminal  is  its  own  obstacle. 
Features  represent  the  fixed  parte  of  the  layout,  the  modules;  wires  represent 
the  flexible  interconnections.  Clearly,  a  sketch  whose  wires  are  writ-behaved 
(e.g.,  consist  of  line  segments)  can  easily  be  encoded  in  a  data  structure. 

2.2.  Legality  and  rout  ability 

I  now  define  what  it  means  to  route  a  sketch,  and  when  a  sketch 
represents  a  legal  layout.  A  curve  in  a  sketch  S  -  (F,  W)  is  a  path  in  the 
plane  that  begins  and  ends  on  features  in  F,  and  intersects  no  features  in 
between.  (Far  example,  the  wires  in  W  are  curve*  in  S.)  Two  curves  in  S 
are  homotopic  if  there  is  a  continuous  deformation  of  the  plane  that  maps 
one  into  the  other  while  holding  the  features  F  fixed.  A  routing  of  S  is  a 


Figure  1.  A  typical  dutch.  Solid  point!  and  Mac  tegmenta  an  ha- 
tore*,  wary  hnea  are  wires,  and  dotted  line*  arc  conceptual  module 
boandariaa. 


aketch  (F,  W')  whose  featurea  are  the  aame,  and  whoae  wirea  can  be  obtained 
from  the  wirea  W  by  a  coutiauoua  deformation  of  the  plane  that  holda  the 
featurea  F  fixed.  (Thua  the  wirea  in  IF'  are  homotopic  to  the  corresponding 
wirea  in  IF.)  A  aketch  ie  aaid  to  be  legal  if  it  repreaenta  a  legal  VLSI  layout, 
which  for  our  purpoaea  meana  that  the  following  condhkma  hold: 

(1)  AH  obatadea  and  wirea  lie  in  the  rectilinear  grid  of  unit  apadng. 

(2)  The  wirea  form  vertex-disjoint  paths  in  the  grid. 

A  aketch  ia  routable  if  it  haa  a  legal  routing.  Uaing  the  algorithm 
in  [6],  a  legal  routing  of  a  routable  aketch  can  be  found  in  polynomial  time. 
Figurea  1  and  2  illustrate  the  concepts  of  legality  and  mutability.  The  aketch 
in  Figure  1  ia  illegal,  because  it  contains  curved  wires.  Nevertheless,  it  is 
routable,  and  one  of  its  legal  routings  is  shows  in  Figure  2. 

2.3.  Routobility  condition* 

The  compaction  algorithm  is  based  on  a  theorem  from  [6]  that  char¬ 
acterises  the  routable  sketches  in  terms  of  the  following  concepts.  If  p  * 
(*p,Vp)  and  q  «  (*,,*«}  are  points  in  the  aketch  S,  then  denotes  the 
open-ended  line  segment  from  p  to  q.  Such  a  segment  is  called  a  cat  if  it 
intersects  no  features  in  5.  The  e opacity  at  a  cut  is  the  w***8"*"*"  number 
of  wires  that  can  legally  crocs  JJ ;  in  symbols, 

cep(H)  «  max{|x,  -  *,J,  |y,  -  y,|,l}  - 1 . 

The  /low  across  P?.  denoted  Jlow{ff),  is  the  number  of  crossings  of  Pf  that 
are  enforced  by  the  topology  of  the  aketch.  (See  Figure  3.)  Crossings  of  pf 
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Figure  1  A  legal  routing  of  the  (ketch  la  Figure  1. 

that  can  be  removed  by  deforming  the  wire*  W  do  not  contribute  to  the  floor. 
More  formally,  ft ov(ff)  is  the  minimum,  over  all  routings  (F,W')  of  (.F, IV), 
of  the  number  of  times  is  crossed  by  wins  in  W. 


Figure  t.  A  portion  of  s  sketch  with  s  cut  ff.  The  flow  across  H  ■ 

I. 

The  notability  of  a  sketch  is  completely  determined  by  the  flows  and 
capacities  of  its  cute.  Let  us  say  that  a  cut  is  taft  [1]  if  its  flow  does  not 
exceed  its  capacity.  Then  we  have  the  following  result. 

Lemma  1.  [6)  A  sketch  that  contains  an  unsafe  cut  is  im¬ 
mutable. 

More  significantly,  the  convene  is  true  (except  when  the  features  of  the 
sketch  are  illegally  placed):  a  sketch  that  contains  no  unsafe  cute  is  routable. 


F.  M.  Matty  *67 

la  fact,  this  statement  may  be  strengthened.  A  critical  cot  ff  is  one  each 
that  p  is  the  endpoint  of  a  feature,  and  q  is  the  closest  point  on  its  feature 
to  p.  The  critical  cuts  are  the  only  important  ones. 

Theorem  3.  [6]  The  sketch  (/*,  W )  is  routable  if  and  only  if 
(F,9)  is  legal  and  every  critical  cut  in  (F,  W)  is  safe. 

The  inequalities  flow(pq)  <  cop[fXj)  far  the  cuts  ??  of  a  sketch  are  called  the 
rout  ability  condition*  for  the  sketch.  Constraints  of  this  sort  will  be  used 
by  the  compaction  algorithm  to  determine  the  optimal  positions  for  layout 
features. 


3.  Computing  flows  in  the  sketch 

This  section  describes  a  procedure  used  to  facilitate  the  computation  of 
mutability  conditions  for  a  sketch.  As  suggested  by  Theorem  2,  the  impor¬ 
tant  attributes  of  a  sketch  are  the  flouts  and  capacities  of  cuts.  Capacities 
are  purely  geometric  quantities,  and  can  be  computed  in  constant  time.  In 
addition,  they  vary  in  a  regular  my  with  the  movement  of  features  during 
compaction.  Flows,  on  the  other  hand,  ate  topological  quantities,  and  are 
relatively  difficult  to  compute.  Moreover,  they  depend  in  complex  ways  on 
the  positions  of  features.  Thus  to  compute  floors,  we  require  a  data  struc¬ 
ture  that  captures  the  topology  of  the  sketch,  and  that  is  invariant  under 
compaction.  I  begin  by  presenting  such  a  structure. 

3.1.  The  adjacency  graph 

The  data  structure  we  need  is  called  the  adjacency  graph  at  the  sketch. 
Its  construction  is  straightforward,  and  is  illustrated  by  Figure  4.  From  a 
point  on  the  leftmost  edge  of  each  obstacle,  except  the  bounding  box,  a  line 
is  drawn  leftward  until  it  hits  another  obstacle.  These  line  segments  and 
rays  will  be  called  hurdle ».  Now  each  wire  is  replaced  by  a  homotopic  wire 
that  intersects  as  few  hurdles  as  possible,  making  sure  that  no  two  wires 
cross.  (Wherever  a  group  of  wires  crosses  a  hurdle  and  crosses  back,  they 
are  moved  to  the  other  side  of  the  hurdle.)  The  resulting  set  of  objects  farms 
a  planar  graph,  whose  vertices  are  obstacles  and  hurdle/wire  crossings,  and 
whose  edges  are  pieces  of  wires  and  hurdles.  The  dual  of  this  graph  is  the 
adjacency  graph  of  the  sketch.  A  node  of  the  adjacency  graph  corresponds 
to  a  face  of  the  original  graph,  and  is  said  to  border  on  the  prints  forming 
the  boundary  of  that  face.  The  adjacency  graph  does  not  change  during 
horisontal  compaction,  because  hurdles  can  only  slide  back  and  forth,  and 
we  will  not  allow  wires  or  features  to  cross  over  one  another. 

The  purpose  of  the  hurdles  is  to  relate  curves  in  the  sketch  to  the 
sketch  topology.  Consider  the  sequence  of  hurdles  crossed  by  a  cum,  fc 
order,  together  with  the  directions  of  crossing.  Such  a  hurdle  eeguent',  can 
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Figure  4.  The  adjacency  graph  of  the  iketeh  in  Figure  1.  Daahed  Hues 
are  hurdles,  and  circles  are  nodes  of  the  adjacency  graph.  Adjacency 
graph  edges  are  not  depicted. 

be  put  into  a  canonical  form  by  removing  all  unnecessary  crossings,  that  is, 
all  places  where  the  curve  crosses  a  hurdle  and  immediately  croeasa  bach  in 
the  other  direction.  One  can  show  that  two  curves  with  the  same  endpoints 
have  the  same  canonical  hurdle  sequence  if  and  only  if  they  are  bomotopic. 

3.2.  Computation  of  flowt 

An  appropriate  search  through  the  adjacency  graph  can  compute  the 
flow  across  a  cut.  To  see  how,  notice  that  there  are  two  kinds  of  edges  in  the 
adjacency  graph:  “wire  edges,"  which  represent  adjacency  across  a  wire,  and 
"hurdle  edges,”  which  represent  adjacency  across  a  hurdle.  A  path  through 
the  adjacency  graph  thus  has  a  hurdle  sequence,  determined  by  the  hurdle 
edges  it  contains.  The  following  lemma  demonstrates  the  correspondence 
between  the  two  kinds  of  hurdle  sequences.  By  the  length  of  a  path  in  the 
adjacency  graph  we  mean  the  number  of  wire  edges  the  path  contains. 

Lemma  3.  Suppose  that  a  cut  of  sketch  5  has  hurdle  se¬ 
quence  H.  Let  Paths  (£?)  be  the  set  of  paths  in  the  adjacency 
graph  of  S,  which  begin  at  a  node  bordering  on  p,  end  at  a  node 
bordering  an  q,  end  have  hurdle  sequence  H.  Then  flow(ff)  is 
equal  to  the  length  of  the  shortest  path  in  Paths  (J5J). 

To  make  use  of  this  lemma,  we  must  be  able  to  find  shortest  paths 
with  given  hurdle  sequences  in  the  adjacency  graph.  It  turns  out  that  the 
following  reedy”  technique  suffices.  We  may  assume  that  each  hurdle  is  to 
be  crossed  from  bottom  to  top.  Call  the  adjacency  graph  G ,  and  define  G' 
to  be  G  with  its  hurdle  edges  removed. 
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Algorithm  F.  (Compute*  the  flow  across  a  cut.) 

Input:  •  cut  with  hurdle  sequence  (hi,...,h«),  the  adjacency  graph  G. 
Local  variables:  integers  t  and  t,  nodes  u  and  v. 

Output:  the  flow  /  across  Jf. 

1.  /  —  min{DlST*FROM(w) :  w  borders  on  p}; 

2.  ftinction  DlST-FROM(w); 

5.  t  —  0;  u  «-  tv; 

4.  for  t  *-  1  to  n  do 
begin 

fi.  v  —  node  bordering  h,-  from  below  that  is  closest  to  «  in  O'; 

6.  t  —  t  +  distance  from  u  to  v  m  O'; 

7.  u  •—  node  adjacent  to  v  across  hurdle  ha 
sad; 

8.  v  «-  node  bordering  q  from  below  that  is  closest  to  u  in  O'; 

8.  return  t  +  distance  from  u  to  v  in  O': 

In  other  words:  for  each  node  bordering  on  p ,  find  the  shortest  path  to  the 
first  hurdle,  cross  the  first  hurdle,  find  the  shortest  path  from  these  to  the 
second  hurdle,  and  so  on.  Breadth-first  search  can  be  used  to  implement 
lines  6-7  and  8-9.  This  approach  may  work  well  in  practice,  but  its  wont 
case  behavior  is  poor.  I  show  in  Section  5  how  to  implement  Algorithm  F 
more  efficiently. 

4-  The  compaction  algorithm 

This  section  defines  mathematically  the  problem  of  compaction  with  auto¬ 
matic  jog  introduction,  and  presents  a  practical  algorithm  that  solves  this 
problem.  I  assume  that  each  layer  of  the  circuit  to  be  compacted  ie  available 
in  the  form  of  a  sketch.  In  fact,  the  compaction  algorithm  will  deal  with  one 
sketch  alone,  for  the  existence  of  multiple  layers  adds  nothing  fundamental 
to  the  compaction  problem. 

4.1.  Configuration  apace 

Let  the  input  sketch  be  denoted  by  5.  For  the  purpose  of  compaction, 
the  obstacles  of  5  should  be  grouped  into  module  §:  collections  of  features 
whose  relative  positions  are  fixed.  The  compactor  is  allowed  to  choose  a 
horizontal  displacement  far  each  module.  Such  a  vector  of  displacements  is 
called  a  configuration  of  S.  The  configuration  d  =  (dj,. . .  ,dn)  corresponds 
to  a  sketch  m  which  module  i  has  been  shifted  right  by  a  distance  d,-  (or 
left  by  a  distance  -d,).  If  the  sketch  has  n  modules,  then  the  set  of  all  its 
configurations  is  the  vector  space  R",  and  the  origin  0  of  this  configura¬ 
tion  apace  corresponds  to  the  original  sketch.  The  sketch  associated  with 
the  configuration  d  will  be  denoted  5(d),  and  is  uniquely  defined  (up  to 
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wire  homotopy)  by  requiring  that  it*  wires  have  the  same  canonical  hurdle 
sequences  they  had  in  the  original  sketch  5(0). 

Using  configurations,  we  can  describe  how  the  sketch  changes  during 
compaction.  If  p  is  a  point  in  5,  its  *  and  y  coordinates  will  be  denoted  zr 
and  yr,  respectively.  The  module  in  which  p  lies  will  be  written  p(p),  so  the 
horisontal  position  of  p  in  the  sketch  5(d)  is  zf  +  d^y  The  notation  p(d) 
stands  for  p  shifted  by  d,  that  is,  the  point  (z,  +  d^j, yr).  In  general,  the 
application  of  an  object  to  a  configuration  d  denotes  the  instantiation  of  that 
object  in  the  sketch  5(d).  Finally,  let  AM(d)  be  difference  in  z-coordinates 
between  q(d)  and  p(d);  that  is, 

AM(d)  =  (zf  +  <*„(,))  -  {*,  +  dM(,j)  . 

4.2.  Problem  ttatement 

The  compaction  problem  is  to  find  a  configuration  d  such  that  5(d) 
is  routable,  and  can  be  routed  in  minimal  width.  As  we  have  stated  it,  the 
compaction  problem  is  generally  very  difficult;  in  fact,  it  is  NP -complete 
[11].  The  reason  is  that  the  routability  conditions  may  not  define  a  convex 
region  of  configuration  space,  and  hence  the  set  of  acceptable  configurations 
(d  €  R"  :  5(d)  is  routable)  is  usually  very  hard  to  search.  For  example, 
whenever  two  modules  are  adjacent,  it  may  be  possible  to  exchange  them 
while  maintaining  the  topology  of  the  sketch,  and  the  resulting  sketch  may 
be  routable.  But  intermediate  positions,  where  the  modules  intersect,  are 
not  routable,  so  the  set  of  acceptable  configurations  is  not  convex.  Even  if  we 
do  not  expect  modules  to  exchange  positions,  loops  of  wire  can  interfere  with 
one  another  in  the  same  way,  as  shown  in  Figure  5.  In  most  optimisation 
problems,  including  compaction,  one  only  expects  to  search  a  convex  sub¬ 
set  of  the  acceptable  configurations,  in  order  to  achieve  a  polynomial-time 
algorithm.  The  algorithm  presented  here  searches  the  largest  such  region 
that  contains  the  initial  configuration,  and  thus  finds  the  best  configuration 
available  to  any  algorithm  of  its  type. 

4.3.  Algorithm  overview 

The  basic  notion  underlying  the  compaction  algorithm  is  that  of  a 
potential  cut.  For  the  purposes  of  this  section,  a  potential  cut  is  a  continuous 
function  that  defines  for  each  configuration  d  €  R"  a  line  segment  between 
two  features  in  5(d).  The  line  segment  may  or  may  not  be  a  cut,  depending 
on  the  positions  of  the  features  in  5(d).  The  configuration  c  is  said  to 
protect  a  potential  cut  i/t  if  either  rl>(c)  is  not  a  cut,  or  ip(e)  is  a  safe  cut. 
The  significance  of  these  definitions  lies  in  a  reformulation  of  Theorem  2  in 
terms  of  potential  cuts.  That  theorem  can  be  read  as  follows: 

Let  5  *  (F,  IF)  be  a  sketch,  and  let  c  be  a  vector  in  its  configura¬ 
tion  space.  There  exist  certain  critical  potential  cuts,  depending 
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Figure  5.  How  wires  can  prevent  modules  from  sliding  past  one 
another.  If  the  upper  module  is  allowed  to  slide  past  the  lower  one, 
the  set  of  acceptable  configurations  is  not  convex. 

only  on  the  features  F,  such  that  the  sketch  5(c)  =  (/'(e),  W (c)) 
is  routable  if  and  only  if: 

(1)  the  sketch  (/’(c),#)  is  legal,  and 

(2)  the  configuration  c  protects  every  critical  potential  cut 
of  5. 

The  compaction  algorithm  works  by  finding  a  subset  of  configuration 
space,  determined  by  simple  linear  inequalities,  whose  configurations  protect 
every  critical  potential  cut.  It  thereby  ensures  that  it  chooses  a  configuration 
that  corresponds  to  a  routable  sketch.  (Condition  (1)  above  can  be  ignored, 
because  the  wireless  sketch  (F(c),#)  is  legal  unless  its  features  fail  to  lie  in 
the  grid.  This  cannot  happen,  because  the  initial  sketch  (F,(V)  is  assumed 
to  be  legal,  and  the  compaction  algorithm  never  considers  nonintegral  dis¬ 
placements  for  modules.)  The  subspace  searched  is  chosen  so  as  to  include 
the  initial  configuration.  An  overview  of  the  compaction  technique  fallows. 

The  central  problem  is  to  find  a  simple  linear  inequality  that  ensures 
that  a  potential  cut,  say  r/>,  is  protected.  One  would  like  to  use  the  routabil- 
ity  condition  eop(^(d))  >  flow  (ip ( d))  as  a  constraint  on  the  configuration  d, 
but  for  most  potential  cuts  V>»  this  constraint  is  not  a  simple  linear  inequality. 
The  difficulty  lies  not  with  the  capacity  of  V>(d),  which  is  determined  solely 
by  the  geometry  of  5(d),  and  depends  in  a  simple  way  on  the  displacements 
d,.  Rather,  the  quantity  flow(rp(d))  is  hard  to  characterise,  because  it  de¬ 
pends  on  the  relation  of  the  line  segment  ^(d)  to  the  topology  of  the  sketch 
5(d). 

The  solution  is  to  find  a  specific  configuration  c  such  that  whenever  the 
potential  cut  v(d)  is  unsafe,  its  flow  is  equal  to  flow(ip(c)).  The  constraint 
eap(ip(d))  >  flow  (ip (c))  is  then  sufficient  to  protect  ip.  Moreover,  when 
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this  constraint  is  written  in  terms  of  the  variables  d,,  it  becomes  a  simple 
linear  inequality,  because  the  right  hand  side  is  constant.  To  find  c,  the 
algorithm  looks  for  a  configuration  that  minimises  the  capacity  of  V>,  subject 
to  the  condition  that  all  critical  cuts  of  smaller  vertical  span  are  safe.  These 
shorter  cuts  force  the  other  features  to  the  side  of  on  which  they  must 
lie,  if  V  is  ever  to  become  unsafe.  If,  in  this  way,  the  algorithm  finds  a 
configuration  c  that  does  not  protect  if),  the  routability  condition  for  ^(c) 
is  remembered.  Otherwise,  the  potential  cut  V>  may  be  ignored. 

4.4.  Deieription  oj  the  compaction  algorithm 

Since  critical  cuts  move  in  nontrivial  ways  during  compaction,  it  turns 
out  to  be  more  convenient  to  consider  the  following  types  of  potential  cuts: 

•  Horisontal  cuts  from  feature  endpoints. 

•  Cuts  between  pairs  of  feature  endpoints. 

The  constraints  generated  from  these  potential  cuts  turn  out  to  be  sufficient 
to  protect  all  the  critical  potential  cuts. 

The  horisontal  potential  cuts  are  particularly  simple,  because  their 
flows  are  independent  of  the  configuration.  These  potential  cuts  are  treated 
first  in  order  to  generate  the  constraints  that  prevent  features  from  crossing 
over  one  another.  A  horisontal  potential  cut  is  a  function  dm  of  the  form 
dM(d)  =  p(d)f  (d),  where  p  and  q  are  points  in  the  original  sketch.  Assuming 
without  loss  of  generality  that  z(  >  xp,  Theorem  2  gives  the  mutability 
condition 

A„(d)  >  /Md*(d))  +  1  . 

Since  flow(tp,[d))  —  flow {pq)  for  pq  horisontal,  the  constraint  is  a  simple 
linear  inequality 

<**»(«)  ”  dMlr)  2  (flo»(P9)  +  1)  +  (*s  -  *«) 

on  the  displacements  of  p  and  q.  The  flow  across  pf  is  easily  computed 
by  Algorithm  F  of  the  previous  section.  The  initial  set  of  constraints  is 
computed  from  the  horisontal  ruts  that  are  incident  on  feature  endpoints. 
Of  course,  the  constraints  are  maintained  as  a  graph  over  the  variables  d,. 
The  constraint  graph  will  be  called  /. 

The  second  stage  of  the  algorithm  concerns  the  cuts  that  are  not  hor¬ 
izontal.  Let  ♦  be  the  set  of  potential  cuts  dm  where  p  and  q  are  feature 
endpoints  with  yp  /  yt .  The  height  of  an  element  dm  of  $  is  the  quantity 
|y,  -  y(|.  This  quantity  is  independent  of  configuration.  Sort  ♦  in  increasing 
order  of  height,  to  form  a  sequence  in  which  flatter  potential  cuts  precede 
taller  ones.  After  dm  €  ♦  has  been  processed,  the  algorithm  can  ensure  that 
the  output  configuration  protects  dm' 
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The  algorithm  examines  the  elements  of  ♦  in  sorted  order,  and  for 
each  one  that  proves  important,  it  adds  an  appropriate  constraint  to  the 
graph  I.  The  constraint  for  a  potential  cut  €  ♦,  with  z(  >  xr,  is  com¬ 
puted  as  follows.  First,  the  algorithm  solves  the  current  constraint  system  I 
together  with  the  temporary  constraint  A„(d)  >  0,  fixing  dM(,),  and  mini- 
mixing  d„(().  Call  the  resulting  configuration  c.  If  c  protects  dH,  then  the 
constraint  set  is  unchanged;  otherwise,  the  constraint 

<**<«)  “  dMlr)  %  (a,  -  *«)  +  flow{+M(c))  +  1 

is  added  to  1.  The  new  constraint  is  a  simple  linear  inequality  derived  from 
the  routability  condition  cop(dM(c))  >  flow(+M(c))  -f  1. 

After  all  the  potential  cuts  in  V  have  been  processed,  the  constraint 
system  I  is  complete,  and  the  algorithm  solves  it  using  a  longeat-path  algo¬ 
rithm.  The  resulting  configuration  is  used  to  build  an  output  sketch,  which 
is  then  routed  using  a  single-layer  router  such  as  Algorithm  R  in  [6].  That 
particular  router  has  the  advantage  of  being  able  to  minimise  the  lengths  of 
the  wires  in  the  touting. 

The  compaction  algorithm  is  summarised  bdow.  We  assume  that  the 
left  and  right  edges  of  the  sketch’s  bounding  box  compose  modules  1  and  n, 
respectively,  and  that  the  top  and  bottom  edges  of  the  box  are  ignored. 

Algorithm  C.  (Compacts  a  sketch  in  the  z-direction.) 

Input:  a  sketch  S  =  (F,  W). 

Local  variables:  the  points  p  and  q,  a  configuration  c,  and  the  constraint 
graph  /  over  variables  d,  (1  <  i  <  n). 

Output:  the  compacted  sketch. 

Subroutines:  Algorithm  F  is  used  to  compute  flows  in  lines  2  and  5;  Dyk- 
stra’s  algorithm  is  used  in  lines  4  and  7;  Algorithm  R  (a  single-layer 
router  horn  [6])  is  used  in  line  8. 

1.  Preprocess  S  as  described  in  Section  3; 

2.  Let  /  be  the  set  of  constraints  {Apffd)  >  flow(pq)  +  1}  where  pf  is  a 

horizontal  cut  with  x,  <  ±%  and  either  p  or  q  is  an  endpoint  of  a 
feature  in  F. 

3.  foreach  pair  of  feature  endpoints  {p,f}  with  xf  <  z(  and  y,  qt  yt,  in 

order  of  increasing  height  do 

begin 

4.  Find  a  configuration  c  that  minimises  cM(t)  -  while  obeying 

the  constraints  1  U  (AH(d)  >  0}; 

8.  If  p(c)q(c)  is  a  cut  in  5(c)  then  if /low(dM(e))  >  cap(d„(c)) 

0.  then  add  the  constraint  AM(d)  >  flow(iM(t))  + 1  to  / 

end; 

7.  Find  a  configuration  c  satisfying  /  that  en  -  e»; 
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8.  Route  the  aketch  5(c)  and  output  the  tenth. 

4.5.  Detaile  of  the  implementation 

•  The  computation  in  line  2  is  easily  performed  using  Algorithm  F  of  the 
previous  section.  The  cuts  themselves  may  be  found  by  any  straightforward 
method,  as  the  algorithm’s  run  time  will  be  dominated  by  other  factors. 

•  Line  3  requires  that  pairs  of  feature  endpoints  be  enumerated  in  order 
of  vertical  separation.  Writing  down  the  pairs  and  sorting  them  would  waste 
large  amounts  of  space;  the  following  approach  is  better.  First  sort  the 
feature  endpoints  by  ^-coordinate,  and  associate  with  each  endpoint  the  next 
higher  endpoint.  Place  these  pairs  in  a  priority  queue,  and  keep  the  queue 
sorted  by  difference  in  ^-coordinates.  At  each  iteration  of  the  loop  (lines  3- 
6),  withdraw  the  best  element  {p,  q }  from  the  priority  queue,  and  process  the 
potential  cut  Then  find  the  next  endpoint  q'  above  q  in  y -coordinate,  if 
one  exists,  and  insert  the  pair  { p,q '}  into  the  priority  queue.  This  method 
uses  liiw  space,  and  no  more  time  than  other  parts  of  Algorithm  C. 

•  To  solve  the  constraint  system  in  line  4,  it  suffices  to  compute  longest 
paths  from  the  vertex  p(p).  Dijkstra’s  algorithm  can  be  used  for  the  purpose, 
because  every  edge  in  the  graph  has  weight  sero  or  less.  (Normally,  Dykstra’s 
algorithm  is  used  to  find  shortest  paths,  and  then  the  edge  weights  must  be 
nonnegative.)  To  see  why  edge  weights  are  nonpoaitive,  consider  the  case 
when  all  the  displacements  d,  are  sero.  Using  the  assumption  that  the  initial 
configuration  is  legal,  one  can  prove  that  it  obeys  all  constraints.  Hence  if 
dj  —  d,  >  ay  is  a  constraint  in  /,  then  it  holds  under  the  assigns  eat  d  =  0. 
The  result  is  that  0  -  0  >  ay,  that  is,  ay  is  nonpoaitive.1 

•  Once  the  algorithm  finds  the  key  configuration  e  in  line  4,  line  5  must 
determine  whether  4>„{c)  is  a  cut,  and  if  so,  what  value  ftow(t„(c))  takes 
on.  The  former  possibility  may  be  checked  by  looking  for  a  feature  that 
intersects  4M(c).  To  compute  /Zow(dM(c)),  Algorithm  F  is  invoked.  It 
requires  as  input  the  hurdle  sequence  of  ^M(c),  which  can  be  found  by 
checking  every  hurdle  that  lies  between  yf  and  yt  in  altitude.  Include  only 


*  This  simple  observation  has  important  implications  for  the  problem  of  com¬ 
paction,  with  or  without  jog  introduction.  Essentially,  it  says  that  whenever  one’s 
input  layout  satisfies  all  the  compaction  constraints,  one  can  ensure  that  all  the 
edge  weights  in  the  constraint  graph  have  the  same  sign,  just  by  choosing  the 
proper  coordinate  system.  The  compaction  variables  must  represent  offttU  from 
the  initial  module  positions,  or  equivalently,  the  reference  point  of  each  module 
must  be  placed  on  the  y-axis.  Then  Dykstra's  algorithm  may  be  used  to  solve 
the  constraint  system,  giving  asymptotic  performance  on  a  graph  G  *  (V,ff)  of 
0(|£|  +  (V  |  log  | V f)  |3),  as  opposed  to  0(|V||E|)  for  the  general  Bellman-Ford  al¬ 
gorithm.  See  (13]  for  a  discussion  of  this  idea. 
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those  hurdles  that  cross  p(e)q(c)  in  5(c).  (A  hurdle  F#  transforms  to  r(c)s(c) 
in  5(c),  just  like  any  other  cut.)  The  hurdle  sequence  should,  of  course, 
be  sorted  by  incoordinate,  and  all  crossings  must  be  from  bottom  to  top. 
Presorting  all  the  hurdles  by  incoordinate  eliminates  the  need  to  sort  each 
individual  hurdle  sequence. 

s  In  line  7,  Dijkstra's  algorithm  should  be  used  once  again,  this  time 
computing  longest  paths  in  /  from  module  1,  which  is  the  left  edge  of  the 
bounding  box  of  the  sketch.  If  desired,  the  designer  or  design  system  may 
add  other  simple  linear  inequalities  to  /,  provided  that  they  are  all  satisfied 
by  the  initial  layout  5(0). 

s  The  configuration  c  found  in  line  7  specifies  the  optimal  compacted 
sketch,  but  that  sketch  must  still  be  constructed  at  line  8.  For  the  purpose 
of  applying  Algorithm  R,  the  single-layer  wire-router  of  [6],  it  is  not  necessary 
to  construct  a  complete  sketch  5(c),  but  only  to  produce  something  called 
the  rvbbtr-band  equivalent  (R.B.E.)  of  5(c).  The  features  of  the  R.B.E.  are 
the  same  as  those  of  5(c),  and  can  be  located  easily.  The  wires  of  the  R.BJ5. 
can  be  found  as  follows.  The  set  of  points  not  lying  on  features  or  hurdles  of 
5(c)  is  a  simply  connected  region,  and  its  boundary  is  polygonal  (if  we  allow 
vertices  at  infinity).  Hence  it  can  be  triangulated  quickly,  and  the  resulting 
set  of  triangles  forms  a  tree  under  the  obvious  adjacency  relation.  We  can 
therefore  find  for  each  wire  tv  in  5(c)  the  shortest  sequence  of  triangles  that 
a  routing  of  iv  could  pass  through,  and  apply  Algorithm  W  from  [6]  to  find 
the  wire  in  the  R.B.E.  corresponding  to  tv. 

4.6.  Complexity  analysis 

We  evaluate  the  time  performance  of  Algorithm  C  in  three  parts.  The 
notation  |D|  is  used  to  mean  the  rise  of  the  data  structure  D.  In  the  worst 
case,  the  most  time-consuming  portions  of  the  algorithm  are  the  0(|F|*)  calls 
to  Dijkstra’s  algorithm  in  line  4.  Each  call  takes  time  0(|£|  +  |V|log  |V|)  on 
a  graph  (V,  £)  [3].  Since  |£|  is  0(|/|),  and  |V|  is  n,  the  number  of  modules, 
this  leads  to  an  estimate  of  0{|F|S(|/|  +  nlogn))  time.  Another  contributor 
to  the  running  time  is  the  call  to  Algorithm  F  in  line  5.  If  Algorithm  F  is 
implemented  as  suggested  in  the  following  section,  then  each  of  its  0(|F|*)  in¬ 
vocations  requires  0(|F|  log  |F|  log  |<?|)  time,  where  G  is  the  adjacency  graph 
at  the  sketch.  So  this  portion  of  the  algorithm  takes  0(|F|*log|F|)og|(?|) 
time.  Finally,  careful  analysis  shows  that  the  construction  and  routing  of 
the  output  sketch  requires  only  0(|F||G|  log  |G|)  time  [6].  Other  operations, 
including  preprocessing  the  sketch,  require  considerably  leas  time  under  the 
reasonable  assumption  that  |W|,  the  number  of  wire  segments  in  the  input, 
isO(|F|*). 

Which  of  the  three  mqjor  parts  of  Algorithm  C  will  dominate  in  prac¬ 
tice  is  not  clear.  In  the  worst  case,  |G|  can  be  as  high  as  n(|F||W|),  if  each  of 
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0(jW|)  wire  segments  internets  a  large  fraction  of  the  bardies,  end  the  cross 
inp  ore  not  redundant.  In  most  situations,  however,  |G|  should  be  closer  to 
|F|.  Making  reasonable  estimates  about  the  average  run  time  of  Algorithm  F 
and  the  density  of  the  constraint  graph  /,  one  can  predict  that  actual  per* 
formance  for  the  entire  operation  will  probably  approach  6(|/’|,+‘)  for  some 
small  value  at  c. 

Space  usage  is  saner  to  evaluate:  the  main  contributors  are  the  graphs 
G  and  /,  along  with  Algorithm  R,  which  may  use  0(|F||£|)  apace  in  the 
worst  case.  Thus  the  worst  case  bound  is  0(|F|*|fF|),  but  none  at  the  data 
structures  of  Algorithm  C  or  Algorithm  R  is  likely  to  approach  its  maximum 
aise.  The  actual  figure  will  depend  on  the  number  of  crossinp  between  wires 
and  **«"  cuts  in  the  sketch  (t.g.,  hurdles),  and  will  probably  be  0(|/’j1+*) 
far  some  constant  o  €  (0,1). 

fi.  Improvements  and  optimisations 

The  compaction  algorithm  has  been  stated  in  its  simplest  farm,  and  there 
is  considerable  room  to  improve  the  performance  of  many  of  its  steps.  This 
section  collects  several  suggestions  for  speeding  up  the  algorithm.  The  mast 
important  of  these  ideas  is  to  add  a  preprocessing  phase  to  Algorithm  C  that 
allows  flows  across  cuts  to  be  computed  quickly. 

5.1.  Fatter  computation  of  flow 

Recall  that  the  compaction  algorithm  computes  flows  using  Algorithm  F 
from  Section  3.  The  most  time-consuming  steps  at  that  algorithm  involve 
searching  through  a  subgraph  G'  of  the  adjacency  graph  of  the  input  sketch. 
Algorithm  F  can  be  implemented  efficiently  by  taking  advantage  of  the  fact 
that  the  graph  G'  is  a  tree.  The  first  task  is  to  pr* process  G'  so  that  one 
can  quickly  determine  the  distance  between  any  pair  of  its  nodes,  and  hence 
speed  up  lines  7  and  9  in  Algorithm  F.  The  second  task  is  to  preprocsss  G' 
so  that  one  can  compute  efficiently  the  closest  node  in  a  connected  subset 
of  G'  to  a  given  node.  This  ability  is  sufficient  to  implement  lines  0  and  8 
of  Algorithm  F,  because  the  set  of  nodes  bordering  a  hurdle  or  feature  from 
below  is  connected  in  G'.  We  wish  to  minimise  the  amount  of  space  need 
by  our  data  structures;  in  particular,  we  cannot  store  explicitly  the  distance 
between  every  pair  at  nodes. 

The  solution  is  to  decompose  recursively  the  graph  G',  and  store  with  a 
node  only  the  distance  to  each  of  its  ancestors  in  the  decomposition.  Figure  6 
shows  the  construction.  Let  n  be  the  number  of  nodes  m  the  graph  G*.  The 
separator  theorem  far  trees  [12]  implies  that  G'  contains  a  verier  r  whose 
removal  disconnects  G'  into  subtrees  containing  st  most  |n  nodes  sach. 
Moreover,  the  vertex  r  can  be  found  in  linear  time  using  depth-first  search 


Figure  6.  A  tree,  drawn  with  solid  Kan,  oad  it*  dwaapodth*  tree. 

to  compote  the  rise*  of  subtree*.  If  we  store  with  every  node  its  distance 
from  r,  then  the  distance  between  two  nodes  m  different  subtrees  of  r  can 
be  computed  by  summing  their  distances  from  r.  Now  we  decompose  the 
subtrees  of  r  recursively,  fanning  a  decomposition  tree  D  an  the  nodes  of  C. 
At  each  stage  of  the  decomposition,  the  rise*  of  subtrees  are  reduced  by  a 
constant  factor,  so  D  has  logarithmic  height.  To  find  the  distance  between 
two  vertices  in  C',  one  first  find*  their  lowest  common  ancestor  (LCA)  in  D, 
sad  then  sum*  their  distances  from  that  ancestor.  This  procedure  clearly 
take*  at  most  O(log  |C|)  time. 

Some  extra  preprocessing  i*  needed  before  we  can  compute  closest 
members  of  connected  seta  at  G.  Let  D  be  the  same  decomposition  tree  as 
above.  The  LCA  in  D  at  a  connected  set  C  £  G'  is  a  member  of  C,  and  we 
can  compute  in  advance  die  LCA’s  of  the  connected  seta  that  we  care  about. 
We  also  store,  for  each  node  y  and  for  sack  of  it*  ancestors  z,  the  highest 
vertex  in  D  that  is  interior  to  the  to  the  path  between  x  and  y  in  G1.  In 
case  z  and  y  are  adjacent  in  G',  we  store  nil  instead. 

I  now  describe  how  to  compute  the  closest  vertex  in  a  connected  set 
C  £  G'  to  a  node  u.  If  u  lies  in  C,  then  the  answer  is  a.  Otherwise,  set 
v  to  be  the  LCA  at  C.  It  v  is  not  an  ancestor  at  u,  then  replace  u  by  the 
LCA  of  u  and  v.  Now  perform  a  "binary  search”  on  the  path  between  «  and 
v  as  follows.  Let  s  be  the  highast  vertex  in  D  between  v  and  *.  It  $  is  nil 
then  the  answer  is  v;  otherwise  setv*-*oru«-»  according  to  whether  s 
does  or  does  not  lie  in  C,  and  repeat  the  process.  After  at  most  0( log  |G|) 
iterations,  the  search  will  terminate,  and  each  iteration  can  be  made  to  run 
in  0(log  |C|)  time. 

5  2.  Compressing  the  eefyeceney  graph 

Both  the  time  and  space  performance  of  Algorithm  C  can  be  improved 
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by  reducing  the  rise  of  the  adjacency  graph.  The  trick  k  to  chooee  hordk* 
in  ouch  »  way  ao  to  rnfakalaa  the  nomber  of  nnooinp  botwai  wire*  aad 
kvdko.  Although  we  defined  hnrdloo  *o  that  way  obotado  kaa  oak  — 
hordk  incident  on  it*  loft,  thia  property  la  unimportant.  Tbo  hordk*  can  bo 
cboaon  to  bo  any  act  of  horiaontal  cate  ooch  that  tba  oat  of  potato  iaaido  tba 
bounding  beat,  bat  not  bring  an  a  hordk  or  a  foature,  k  limply  connirtad 
Efahakntbi  it  obetack*  and  hordko  arc  coaoidorad  aa  tbo  no  doe  aad  adgaa, 
respectively,  of  a  graph,  then  this  graph  mot  bo  a  tree. 

We  can  aao  a  minim-coat  apannmg  tree  algorithm  to  Ind  a  oat  of 
hordko  that  crooa  aa  few  vino  ao  poaoibk.  Every  boriaontal  cat  botaroon 
different  obotacka  k  a  potential  hordk,  bat  we  may  reotrict  oar  attention 
to  boriaontal  cate  that  are  incident  an  baton  endpointo.  There  are  at  moot 
0(|-F|)  each  cote,  and  they  con  be  thought  of  oo  the  odgao  of  a  graph  B  over 
the  obstacles.  The  coot  of  an  edge  will  be  the  nambor  of  nnaaiwp  of  the 
cat  by  wireo  in  the  original  aketch;  coata  can  be  coapated  effkkntly  aaing 
a  ac arming  algorithm  aa  in  (6j.  The  hordko  are  cboaon  to  bo  the  odgaa  in  a 
minimum-coat  ■p»««*««g  true  of  the  graph  B. 

5.3.  Ignoring  unimportant  potential  eats 

A  few  aimpk  chocka  will  probably  apeed  op  Algorithm  C  by  a  factor 
of  4  or  more.  Firat  of  all,  a  corollary  of  the  corroctneaa  proof  k  that  a  po¬ 
tential  cot  whooe  charity  k  minimal  in  the  initial  configuration  cannot 
generate  a  cautraint.  Therefore,  the  algorithm  need  only  check  potential 
rate  4m  for  which  |sf  -  *p\  >  (yf  -  yp\.  Second,  the  loww  endpoint  of  a 
feature  need  not  be  conoidered  in  conjunction  with  feature  endpointo  above 
it,  and  symmetrically  for  the  upper  endpoint.  Similarly,  a  potential  cut 
which  in  the  initial  configuration  travak  rightward  from  the  kft  endpoint 
of  a  boriaontal  feature  need  not  be  conoidered,  aad  symmetrically  for  right 
endpointo.  Finally,  a  potontial  cut  with  st  >  a,  need  not  be  chocked  if 
in  all  configuration  d  each  that  AH(d)  >  0,  the  line  aegment  dn(d)  k  not 
a  cut.  Thia  k  known  to  be  true,  for  exampk,  if  every  rightward  cut  Cram  p 
croeeee  a  feature  in  the  aame  moduk  ao  p. 

0.  Highligbto  of  the  corroctneaa  proof 

The  correctneoo  proof  of  Algorithm  C  cannot  be  included  here  due  to  apace 
reetrictiona,  for  k  k  highly  involved.  It  does,  however,  lead  to  aa  uader- 
atandmg  of  the  aeeumptiona  that  underlie  the  compaction  algorithm.  In  thk 
•action  I  etate  a  aet  of  condition  under  which  one  can  contract  a  com¬ 
paction  algorithm  that  perform*  automatic  Jog  introduction  in  the  mittnir 
of  Algorithm  C.  Knowing  theee  condition,  one  can  determine  whether  Al¬ 
gorithm  C  will  continue  to  work  if  the  definition  of  a  legal  aketch  k  modified. 
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We  require  t  marc  technical  definition  of  potential  ait.  Let  P  md 
Q  be  features  in  the  original  sketch  5.  and  let  9  be  a  continuous  fraction 
that  defines,  for  each  configuration  d,  a  Mm  eegmeat  between  the  features 
P(d)  and  0(d)  fa  the  (ketch  5(d).  The  fraction  9  fa  a  potential  cat  if  tha 
position  o {  9(d)  relative  to  P(d)  and  9(d)  depends  only  on  tbe  djcplac— cat 
between  P(d)  and  9(d),  aaadjr  A/>g(d)  ■  4^q)  -  d,(p)  (Btretebiag  tha 
notation  slightly).  In  other  word*,  9  moot  satisfy  the  following  condition. 

If  d  and  d1  are  any  two  configuration*  *och  that  Apg(d)  * 
Aro(d'),  then  9(d')  i*  equal  to  9(d)  ahifted  to  tbe  ri^ht  by 
*,(r)  ~  d*r)  *aiu- 

Algorithm  C  compact*  a  aketch  5  by  computing  a  sequence  #(S)  ■ 
(9  i,...,9«.)  of  potential  rate  of  S,  and  wring  their  routability  condition*  to 
constrain  the  position*  of  module*  in  5.  Tbe  proof  that  thi*  technique  fen* 
crate*  a  sufficient  and  optimal  constraint  system  depends  on  several  facts. 
One  each  fact  is  Lemma  1,  which  ensure*  that  an  unsafe  cut  makes  a  sketch 
unrou table.  Another  important  fact  is  that  the  potential  cots  9(5)  deter¬ 
mine  the  routability  of  the  modified  sketches  5(d).  Specifically,  they  have 
the  following  property: 

Rovtokiktg  property.  If  5(0)  is  routable,  and  lot  aB  A  €  [0,1]  the 
configuration  Ad  protects  every  9  €  0(5),  then  5(d)  is  routable. 

The  capacities  of  the  potential  cuts  also  have  a  special  property: 

Bitonic  property.  For  each  d  €  t(f),  and  each  fine  L  in  con¬ 
figuration  space,  there  is  a  point  c  of  L  at  which  the  capacity 
cap(9(c))  is  minimal,  and  cap(9(d))  is  nonderroasnig  as  d  moves 
away  from  c  along  L. 

la  principle,  our  compaction  method  depends  on  only  one  further  fact. 


Ordering  property.  Suppose  that  the  following  statements  hold. 

(1)  The  configuration  d  protects  9<  for  all  •  <  k. 

(2)  The  configuration  d  lies  on  tbe  boundary  of  the  set  (c  € 

R*  :  9(c)  is  a  cut). 

(3)  The  cut  s  is  properly  contained  in  the  line  — g"»—«  9*(d). 

Then  s  is  safe  in  the  aketch  5(d). 

Algorithm  C  satisfies  this  requirement  by  sorting  its  potential  cuts  according 
to  height. 

The  conditions  listed  above  allow  us  to  determine  when  Algorithm  C 
can  be  extended  to  more  general  kinds  of  layouts.  Suppose  that  under  new 
definitions  of  "■ketch”  and  "legal,”  Lemma  1  continues  to  hold,  and  that 
given  a  legal  sketch  5,  one  can  compute  in  polynomial  time  a  sequence 
f  (5)  of  potential  cuts  with  the  routability,  bitonic,  and  ordering  properties. 


Suppose  also  that  ooe  hai  a  angle-layer  router  that  can  £11  ia  the  vine  of  5 
after  compaction.  Then  the  technique  c f  Algorithm  C  gives  a  polynomial¬ 
time  algorithm  to  compact  the  aketch  S. 

7.  Extensions  and  discussion 

This  section  suggests  several  ways  in  which  the  compaction  algorithm  can 
be  extended,  and  concludes  with  an  overview  of  the  main  results. 

7.1.  Generalisations  of  Algorithm  C 

7.1.1.  Multiple  layers.  We  have  assumed  that  the  layout  to  be  compacted 
comprises  only  one  layer.  To  remove  this  restriction,  Algorithm  C  need 
only  generate  constraints  for  each  layer  of  the  design  independently,  as  in 
lines  1-6.  Since  most  modules  have  components  on  more  than  one  layer,  the 
resulting  constraint  systems  must  be  merged,  by  choosing  the  moat  restric¬ 
tive  constraint  between  every  pair  of  modules.  The  merged  system  is  then 
solved  normally. 

7.1.2.  Wire  widths.  None  of  my  results  depend  on  the  assumption  that 
all  wires  have  the  same  width.  If  they  do  not,  one  need  only  replace  the 
counting  of  wires  by  the  summation  of  their  widths,  in  the  definition  of  flow, 
in  Lemma  3,  and  in  Algorithm  F.  With  suitable  modifications,  features  can 
have  differing  thicknesses  as  well.  More  details  are  included  in  [14]. 

7.1.3.  Multiterminal  nets.  One  major  deficiency  of  the  aketch  as  a  layout 
model  its  limitation  on  the  number  of  terminals  that  may  contact  a  wire.  The 
model  can  be  generalised,  however,  and  the  compaction  algorithm  remains 
essentially  unchanged.  (One  difference  is  that  Algorithm  R,  the  single-layer 
router,  can  no  longer  minimise  the  wire  lengths  in  the  output  aketch.)  See  [6] 
for  a  discussion  of  how  to  handle  multiterminal  wires. 

7.2.  Summary  and  conclusion 

The  main  theoretical  contribution  of  this  paper  is  a  polynomial-time 
algorithm  that  compacts  IC  (or  PCB)  layouts  while  introducing  jogs  into 
wires  in  an  optimal  fashion.  The  power  of  Algorithm  C  comes  from  the 
elimination  of  wires  as  hard  objects  in  the  layout,  and  their  replacement 
by  constraints  between  modules.  The  use  of  routability  conditions  to  solve 
placement  problems  is  not  new  [7,15,16],  but  until  now,  only  channel  routing 
problems  had  been  considered.  The  reason  is  that  the  mutability  of  general 
planar  layouts  was  not  adequately  understood  until  very  recently  [1,6].  To 
characterise  planar  mutability  requires  a  robust  model  of  a  circuit  layer, 
such  as  the  aketch,  and  a  fair  amount  of  theory.  In  addition,  some  care  is 
needed  to  apply  mutability  conditions  to  the  compaction  of  general  sketches; 
proving  the  correctness  of  Algorithm  C  is  nontrivial. 
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On  the  practical  aide,  my  compaction  method  can  be  expected  to  pro¬ 
duce  high-quality  layout*  with  little  designer  intervention,  saving  both  in 
chip  area  and  design  time.  Its  primary  drawback  lies  in  its  nse  of  com¬ 
putational  resources.  Although  there  are  good  reasons  to  believe  that  its 
wont-case  performance  bounds  will  not  be  approached  in  practice,  resource 
limitations  may  prevent  it  from  being  used  to  compact  large  layouts  all  at 
once.  Algorithm  C  is  amenable  to  use  at  all  levels  of  the  design,  however,  so 
that  hierarchical  compaction  can  alleviate  much  of  the  resource  problem.  It 
also  may  be  suited  to  use  in  channel  routing,  where  the  number  of  compo¬ 
nents  is  not  too  great.  The  idea,  which  was  implemented  at  BeU  Labs  (see 
Acknowledgements)  is  as  follows:  the  channel  is  artificially  inflated,  so  that 
an  an  ordinary  channel  routing  algorithm,  which  may  have  difficulty  with 
crowded  channels,  may  succeed;  then  a  compactor  like  Algorithm  C,  with 
the  ability  to  insert  arbitrarily  complex  jogs,  is  applied  in  order  to  compact 
the  channel  back  to  the  proper  aise. 

One  important  question  left  open  by  my  research  is  whether  the  com¬ 
paction  method  embodied  in  Algorithm  C  is  more  efficient  in  practice  than 
the  straightforward  algorithm,  namely,  inserting  jog  points  into  each  wire 
where  it  crosses  each  horisontal  gridline,  and  solving  the  resulting  constraint 
system  normally.  This  technique  is  evidently  simpler  than  that  of  Algo¬ 
rithm  C,  and  may  be  man  efficient  in  practice.  On  the  other  hand,  it  may 
be  possible  to  extend  Algorithm  C  to  situations  where  wires  and  modules 
may  contain  diagonal  segments,  and  grid-based  algorithms  break  down. 
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